Line-broadening in gravitational radiation from gamma-ray burst-supernovae 



O 

o 



Maurice H.P.M. van Putten 
LIGO Project, NW17-161, 175 Albany Street, Cambridge, MA 02139 

Hyun Kyu Lee and Chul H. Lee 

Department of Physics, Hanyang University 133-791, Seoul Korea 

Hongsu Kim 

Astronomy Program, SEES, Seoul National University, Seoul 151-742, Korea 

Core-collapse in massive stars is believed to produce a rapidly spinning black hole surrounded by 
a compact disk or torus. This forms an energetic MeV-nucleus inside a remnant He-core, powered by 
black hole-spin energy. The output produces a GRB-supernova, while most of the energy released 
is in gravitational radiation. The intrinsic gravitational-wave spectrum is determined by multipole 
mass-moments in the torus. Quadrupole gravitational radiation is radiated at about twice the 
Keplerian frequency of the torus, which is non-axisymmetric when sufficiently slender, representing 
a "black hole-blob" binary system or a "blob-blob" binary bound to the central black hole. We here 
discuss line-broadening in the observed spectrum due to Lense-Thirring precession, which modulates 
the orientation of the torus to the line-of-sight. This spectral feature is long-lived, due to weak 
damping of precessional motion. These events are believed to occur perhaps once per year within a 
t-H ' distance of lOOMpc, which provides a candidate source for Advanced LIGO. 

m : 

y—i . LIGO recently completed its second science run 0. The results demonstrate the feasibility of km-sized laser 
interferometric gravitational-wave detectors. In the broad band width of 20-1000Hz, advanced LIGO is sensitive to 
extragalactic events of binary neutron-star coalescence and GRB-supernovae from rotating black holes. This band 
width also covers emissions from new born neutron stars, binary coalescence of black holes, as well as emissions 

^vq i from rapidly rotating neutron stars. Anticipating the spectrum of gravitational radiation from candidate sources is 

important in designing optimal detection strategies. While the spectrum of binary inspiral is well-understood to high 

post-Newtonian order at large binary separations, that of GRB-supernovae is only beginning to be identified. 

Core-collapse of massive stars is considered a potential site of gravitational radiation 0,0, H, @ • Here, 

gravitational radiation is produced by the release of gravitational binding energy during collapse and in accretion 

Oh' processes on a newly formed black hole (e.g, ||). If accretion is driven by angular momentum loss in magnetic 

winds, the timescales in core-collapse by magnetic regulated hyperaccretion are relatively short: about 0.05s for the 

outer accretion disk, and about 0.01s for the inner accretion disk [lflj . Hyperaccretion flows are probably strongly 

turbulent, which would imply a broad gravitational- wave spectrum. Aforementioned studies appear to indicate an 

energy output, which leaves a range of detectability by ground based detectors of up to about lOMpc. These events 

. should therefore be considered in the context of core-collapse events independent of the GRB phenomen on, in light 

of current estimates on the local GRB event-rate of about one per year within a distance of lOOMpc 1 1 lL Il2|. These 

studies on gravitational radiation in core-collapse of a massive star do not invoke the spin-energy of a newly formed 

' black hole. 
_ _ i 

The lifetime of rapid spin of the black hole surrounded by a torus in suspended accretion is about tens of seconds 
[lf]lll3| . This represents the timescale of "unseen" dissipation of black hole-spin energy in the event horizon, the rate 
of which is governed by the poloidal magnetic field-energy in the surrounding torus 10, 13]. This timescale is in good 
agreement with the durations of long GRBs, which reflect the lifetime of activity of the inner engine [lil ]. 

Optimal detection strategies for gravitational radiation from GRBs may be facilitated by a priori knowledge from 
a specific model. Our model for GRB-supernovae from rotating black holes describes the formation of a GRB 
accompanied by a radiatively-driven supernova, while most of the energy output is in gravitational radiation [l3j . 
The lowest order quadrupole moment is manifest as a "black hole-blob" binary system or a "blob-blob" binary bound 
to the central black hole when the torus is sufficiently slender, giving rise to |L3] 

<w««*«M&)(^)- {mT (1) 

in gravitational radiation. Here, r\ = £It/^h denotes the ratio of the angular velocity fix of the torus to the angular 
velocity of a rapidly rotating black hole of mass Mh ■ This energy output results from a causal spin-connection 
to the black hole of a surrounding torus with multipole mass-moments. For an intrinsic mass-inhomogeneity SMt in 
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the torus, we have a luminosity of gravitational radiation according to 

L gw = {32/5)(M H /R) 5 (SM T /M H f, (2) 

where M. ~ Mh{SMt/Mh) 3 ^ 5 denotes the chirp mass. The relative mass-inhomogeneity 5Mt/Mh is determined 
self-consistently by the balance equations of energy and angular momentum flux in a suspended accretion state to be 
about a few promille \\?\ . The fractions of black hole-spin energy radiated in various channels is independent of the 
mass of the torus (provided Mt is a few tenths of a solar mass). The mass of the torus only affects the lifetime of 
rapid spin of the black hole and hence the duration of the burst. 

The energy in gravitational waves Q is larger than the true energy E 1 in GRB-afterglow emissions (see further 
0), i.e.: 

E gw — 10 3 i? 7 . (3) 

We attribute the GRB-afterglow emissions to baryon poor outflows along an open magnetic flux-tube subtended at a 
finite opening angle on the event horizon of the black hole. Statistical analysis of GRBs with individually measured 
redshifts points towards strongly beamed GRB-emissions [Jjj in the form of a highly anisotropic beam accompanied 
by very weak GRB-emissions over essentially all angles [l2j | . This is distinct from purely conical outflows with uniform 
emissions over a finite opening angle and zero emissions outside. The non-uniform strongly anisotropic beam gives 
rise to an anti-correlation between inferred opening angle and redshift, given a standard true GRB-energy output 
in response to a finite flux-limit of a detector. The weak "all-angle" emissions are proposed to be at the level of 
GRB980425, or about 10 -4 times the standard GRB luminosities - weak, but significant in allowing nearby GRBs to 
be identified off-axis as apparent anomalous, low-luminosity events. 

A radiatively-driven supernova energy accompanies the GRB and burst in gravitational waves powered by magnetic 
winds coming off the torus in suspended accretion around the central black hole |13| . The predicted kinetic energy 

E SN 2 x 10 51 erg (4) 

is in good agreement with the observed energy of about 2x10 ergin the aspherical explosion of SN1998bw 0. We 
note that the "hypernova" energies ^3 °f a f ew times 10 52 erg |l8l are isotropic equivalent energies assuming 
spherical symmetry, not true kinetic energies in the aspherical explosions. 

We here study the spectral signature of Lense-Thirring precession in Q . Lense-Thirring precession can be described 
by nodal precession [20j. In the slow-rotation limit, its angular velocity £Ilt (2Z| agrees essentially to the frame- 
dragging angular velocity. The orientation of the torus is well-defined, in that misalignment between its angular 
momentum and angular velocity is of second order, 0{9VLlt) 1 where denotes the wobbling angle. This contrasts 
with rigid body motion, where precession is due to a misalignment between the angular momentum and angular 
velocity vectors. Lense-Thirring precession is frequently invoked in modelling QPOs produced by accretion disks in 
X-ray binaries j2l|. It represents a general relativistic effect due to frame-dragging, and hence it acts universally on the 
large-scale structure of a torus surrounding a black hole. A torus which is misaligned with the spin-axis of the black 
hole precesses with essentially the frame-dragging angular velocity described by the Kerr metric. In Boyer-Lindquist 
coordinates, we have FIlt — 2Jh/R 3 for a black hole angular momentum Jh = M 2 sin A in terms of the mass Mu 
and the specific angular momentum sin A = a/M. Given the angular velocity £Ih = tan(A/2)/2M of the black hole 
and the angular velocity fix — M 1 / 2 ./? -3 / 2 of the torus, we have 

SiI,2xlO-^) 2 ^(A/ 2 ). (5) 

At nominal values rj ~ 0.1, SIlt is about 10% of FIt, or, equivalently, about 1% of fijy. The associated precession 
of the black hole, due to conservation of total angular momentum, can generally be neglected when the angular 
momentum of the torus is much less than that of the black hole. 

Our study considers Lense-Thirring precession in response to an initial misalignment of the torus relative to the 
spin-axis of the black hole. The active nucleus is believed to form in core-collapse of a massive stars |22) inside a 
remnant envelope (e.g. |l3j|). The angular momentum of the black hole is commonly attributed to orbital angular 
momentum, received during spin- up in a common envelope phase with a companion star [TtI |23| . A misalignment 
in the nucleus may hereby reflect misalignment between the spin of the progenitor star and the orbital motion of 
the companion star. Lense-Thirring precession modulates the orientation of the torus to the line-of-sight, introducing 
phase-modulation in the line-emissions as observed at the detector. This broadens the emission lines, and may produce 
side-bands similar to those arising in amplitude modulation. There arises an additional pair of low similar to radiation 
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coming off a freeiy processing neutron star. These low-frequency emissions produce relaxation, and gradually bring 
the torus into alignment. Their low-luminosities allow the misalignment and its precession to be long-lived. 
For the two polarizations of gravitational waves, we have 

h+ = 2(1 + cos 2 l) cos(2il T t), h x = -4 cos t sin(2f2 T i) (6) 

where i denotes the angle between the angular momentum and the the line-of-sight. Here, we have suppressed a 
coefficient 2 (et/ Mh) [Mh/t) , where ek denotes the kinetic energy of the mass-inhomogeneity m in the torus, and r 
denotes to the distance to the source. Precession of the torus introduces a time- varying angle i(t), which modulates 
the strain amplitudes h+ and h- at the observer (Fig. 1). Given a mean angle to of the angular momentum of the 
torus to the line-of-sight and a wobbling angle 9, the time-dependent angle i(t) of the same satisfies 

cos = sin to sin 8 cos(flLTt) + cos lq cos 8. (7) 

Substitution of cos i(t) into © produces phase-modulation at the Lense-Thirring frequency. The expressions © 
hereby produce line- broadening by phase- modulation about the carrier frequency 2Qt, according to 

h + = 6 cos(2ft T t) + cos(2Q T i + 2t) + cos(2Q T t - 2i), (8) 

and 

h x = -2 [sin(2fM + t) + sm(2Q T t - i)] . (9) 

The resulting wave-forms are illustrated in Fig. 2. 

In case of 8 « i , we may linearize © h+ = h { ° ] + Oh^ + 0(8 2 ), h x = /i ( x 0) + 8h { ^ + 0{8 2 ), 
where h+ = h+(io) and h^p = h x (io) refer to the strain-amplitudes © with t = Lq, and h+ = 
-sin(2<.o)[cos([2fi T + n iT ]t) + cos([2fi T -O iT ]t)], and h^ = 2sin(t ) [sin([2fi T + tt LT ]t) + sin([2fi T - fl LT }t)} . 
For small wobbling angles, therefore, phase-modulation reduces to amplitude modulation, producing a carrier fre- 
quency accompanied by two side-bands at frequencies 

2n T -n LT , 2n T , 2VL T + n LT . (10) 

The ratio of the amplitude line-strengths of the side bands to that of the carrier satisfies 

1/2 

ts of 1 + cos 2 i \ ■ d-n 

K~6 — - 2 ■ — smt , (11) 

\ 1 + 6 cos^ io + cos 4 l J 

where we used CIlt << 2ilr- Averaged over all angles to, we have K ~ 8/2. A wobbling angle of about 30° typically 
produces side-bands of relative strength 20% (taking together h+ and hx in each side-band). 

Low-frequency emission lines are found by considering the projections of a processing torus onto the celestial sphere, 
as shown in Fig. 1. Translating the geometrical picture of [32j to the case of a torus around a black hole, we observe 
that (1) projections along the spin-axis of the black hole on the celestial sphere produces an ellipse which rotates at 
twice the Lense-Thirring frequency; (2) projections of the torus along directions in the orbital plane on the celestial 
sphere produce a line-segment which oscillates at the Lense-Thirring frequency. A precessing torus hereby introduces 
a spectral anisotropy, whereby emissions at twice the precession frequency are preferentially along the spin-axis of the 
black hole, and emissions at the precession frequency are preferentially along directions in the equatorial plane. 

The strain-amplitude of the low-frequency emissions satisfies j3^ 

h LT = sin(20)%^h LT , (12) 



where 



and 



2r 



hLT = e+ ( h\ tan 6 + h\\ + e x (h\ tan 8 + h x ) (13) 



h l l = 2(1 + cos 2 i)cos{2tt LT t),h l l = 4 cos i sin(2fiz, T *), (14) 



h\ = -sin(2t)cos(OiTi))^x = sin tsin(nLT^)- (15) 
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The amplitude ratio of lll.'jll relative to © satisfies 



2 1 

h LT 1 / Q. LT \ I m 



h 8 V fir 7 V U 



•it 



< 1%. (16) 



Here, we note that Mx/m is a ratio of order unity, for m on the order of a few permille of Mjj in a torus of a few 
tenths of M©. Emissions by (|13|l along the axis are second order in the wobbling angle 9; emissions in directions along 
the orbital plane are first order in 9. This reflects the amplitude of the time-variability of the quadrupole moment in 
the association projections on the celestial sphere. The associated net luminosity satisfies 

L £ = ^^t(A/) 2 sin 2 (20)(l + 16 tan 2 9), (17) 

where the first and second part in the bracket correspond to the emissions by h x and h", respectively at frequencies 
£Ilt and 2Q.lt- Their contributions are equal when tan 9 = 1/4 or 6 ~ 14°, and the emissions in h dominate for 
wobbling angles 9 < 14°. In this event, the spectrum produced by (|13fl varies according to the different dependencies 
on the viewing angle t. Viewing the source edge on (the observer in the orbital plane) sees h)\ at 2£1lt and hj, at 
Q-lt- Observers along the spin-axis of the central black hole see no lines at ^lt m ^+ an d h^, but a single line at 
2SIlt in both polarizations of /i" which are second-order in 9. 

The lifetime of precession can be calculated by considering the backreaction of gravitational radiation on the 
torus. Gravitational radiation backreaction forces consist of dynamical self-interactions and radiation-reaction forces 
[24I Eil l3l| . For sources with nonrelativistic motions in the approximation of weak internal gravity, the latter can be 
modelled by the Burke-Thorne potential as it arises in the 2^ post-Newtonian approximation. There is no change 
in the continuity equation in this intermediate order (such as in the second-order post-Newtonian approximation 

HHHIH). 

Gravitational radiation backreaction on an intrinsic quadrupole moment produced by the m = 2 Papaloizou-Pringle 
instability, is governed by the small parameter j3 = (7ri?E/10) (2HtR) 5 , or 



0.1/ V°- 1M o/ V b ) \ m h 

where b ~ M denotes the minor radius of the torus in suspended accretion. This backreaction stimulates the m = 2 
instability [sol ]. 

Gravitational radiation backreaction to the precessional motion is due to the low-frequency emissions, since the 
emissions due to intrinsic multipole mass-moments leave the wobbling angle unaffected. The backreaction of the 
low-frequency emissions are similar but not identical to that on a freely precessing neutron star. Adapting |25| to 
the present case, the low-frequency gravitational radiation reaction torque Tjj£ is perpendicular to the symmetry 
axis of the torus, along which there exists a discrete symmetry in the presence of intrinsic multipole mass-moments. 
Gravitational radiation reaction torques associated with the latter act along this symmetry axis. It extracts the 
angular momentum along the spin-axis of the black hole - the axis about which the torus wobbles - and angular 
momentum associated with the wobble angle itself. In the suspended accretion state, the angular momentum along 
the spin-axis of the black hole is constant on timescales less than the lifetime of rapid spin of the black hole - the 
GRB timescale of tens of seconds. We may follow [2^, by imposing the constraint 

J z = const. (19) 

Reduction of the wobbling angle reflects extraction of the kinetic energy of the wobbling motion of the torus. According 
to j2{|, we have 

T^cos9 

9 = ^-j , (20) 

where T^Jr^i^sin^ = [2q. At late times, wobble angle is small, whereby 9 = —9/tq with Tg 1 = 

(1/10) {VLlt/^hY £1% Mt R 2 for a torus of mass Mt and major radius R. In dimensional units, and using JSJ, 
we have 

r] \-28/3 / M H 
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For a wide range of values of r\ and Mh, the damping time due to the radiation backreaction is longer than the 
time-scale of tens of seconds of long GRBs. 

We may compare the strength of the high and low-frequency emissions. The luminosity at frequency 2Qj> due to 
an intrinsic quadrupolc moment in the torus of mass My satisfies 

■ rj \ 10 / 3 ( m y 



L gw (2Q T ) = 3 x 10- (^j 1 — 1 • (22) 



The luminosity in emissions of h associated with a small wobbling angle satisfies 

L gw (Q LT ) = 1.6 x 10~ 9 (^) 10/3 (H) V (23) 

where we used Qlt/^t — V- ln the suspended accretion state, we expect a lumpiness m ~ 0.15%Mg or larger, 
whereby 

It follows that 



2 

0.1 1 \0.1M T 



L gw (n LT ) : L gw (2n T ) = 5 x 10-» [ — ) I I . (25) 



We conclude that the low-frequency emissions are low- luminosity lines, relative to the expected emissions from an 
intrinsic quadrupole moment. This low-luminosity is due to the strong frequency dependence of gravitational radiation, 
and the limited source of (kinetic) energy in the wobbling motion of the torus. In contrast, the gravitational radiation 
produced by the intrinsic quadrupole moment is produced by the Keplcrian frequency and it derives from the large 
reservoir in spin-energy of the the black hole. 

To summarize, we predict long bursts of gravitational radiation in GRB-supernovae from rotating black holes with 
an energy output of 4 x 10 53 erg. These bursts represent the catalytic conversion of black hole-spin energy into gamma- 
rays and gravitational waves, and an accompanying supernova. The gravitational wave-luminosity is dominated by 
emissions from multipole mass-moments in the torus for sufficiently slender tori (minor-to-major radius of the torus 
less than about 0.75 for a mass-moment m = 1, and less than about 0.33 for mass-moment m = The duration 

is set by the lifetime of rapid spin of the black hole. A misalignment of the torus and the black hole conceivably 
represents a remnant feature of the angular momentum distribution in the progenitor binary of a massive star and 
its companion. It produces phase-modulation in the observed strain-amplitude. 

A detailed understanding of the wave-form of gravitational radiation is important in applying matched filtering. 
For matched filtering, we have the signal-to-noise estimate 

s\ / g^ooHz; i y 1 3 /2 5/2 

NJ mf \4x lO-^Hz" 1 / 2 / 101 78 K ' 

Presently, the second science run in LIGO shows s£ /2 (500Hz) = 5 x lO^Hz- 1 / 2 I], which corresponds to S/N = 8 
for a source at ~ IMpc. This will improve to a sensitivity range of lOOMpc with the Advanced LIGO noise-strain 
amplitude target Si ' 2 (500Hz) = 4 x 10~ 24 Hz~ 1,/2 . This sensitivity range corresponds a true event-rate (seen in 
gravitational radiation; seen and unseen in GRB-cmissions) of one per year, based on a GRB beaming factor of about 

500 GHE2- 
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FIG. 1: The geometry of Lense-Thirring precessing torus with angular velocity 51t around a rotating black hole with angular 
velocity Qh- The wobbling angle is 6. The line-of-sight has an angle to relative to the spin-axis of the black hole. The angle of 
the orientation of the torus to the line-of-sight varies in time, between to — 9 and to + 0. This introduces phase-modulation in 
the observed strain-amplitude of emissions by an intrinsic quadrupole moment in the torus. Low-frequency spectral anisotropy 
is produced by apparent quadrupole moments, corresponding to projections 1T\\ and %± of the torus onto the celestial sphere in 
directions along and orthogonal to the spin-axis of the black hole. The associated modulations are at twice and, respecively, 
once the Lense-Thirring precession frequency. 
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FIG. 2: The strain-amplitudes h+ and h x around twice the Keplerian angular frequency fir are produced by a quadrupole 
mass-moment in the torus in a state of suspended accretion. The observed strain-amplitudes are subject to phase-modulation 
by Lense-Thirring precession of the torus. For wobbling angles which are small relative to the mean angle to the line-of-sight, 
the resulting amplitude modulation produces two side-bands to the carrier frequency 2Qt, separated by the Lense-Thirring 
frequency Qlt- Shown are the strain-amplitudes for the case of to = 0, 7r/8, 7r/4, tt/2 for a wobbling angle of 30° and a Lense- 
Thirring precession frequency of 1/8 of the Keplerian frequency. The amplitude corresponds to a source at unit distance, and 
the index refers to the number of Keplerian periods. Additionally, small contributions to the strain amplitude arise at once 
and twice Qlt in response to apparent quadrupole moments in the projections of the torus onto the celestial sphere. 



